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We discuss consequences of Fermi exchange symmetry on collisions of polar molecules at low 
temperatures (below 1 K), considering the OD radical as a prototype. At low fields and low tem- 
peratures, Fermi statistics can stabilize a gas of OD molecules against state-changing collisions. We 
find, however, that this stability does not extend to temperatures high enough to assist with evap- 
orative cooling. In addition, we establish that a novel "field-linked" resonance state of OD dimers 
exists, in analogy with the similar states predicted for bosonic OH. 



I. INTRODUCTION 

Cooling and trapping molecules in their ro- vibrational 
ground states has proven to be a daunting experimental 
task, yet now it has been achieved Q, H0]. Currently, the 
samples produced by Stark slowing are limited to tem- 
peratures around 1 mK, which is cold enough to trap, 
but not yet cold enough for interesting applications to 
novel dilute quan tum gases of fermions [E IE IE IE El 
or bosons |E HE El El EE 13 To produce colder, 
denser samples for these applications, an attractive ap- 
proach may be to use sympathetic cooling with an eas- 
ily cooled species (such as Rb) 15], or else evaporative 
cooling. Knowledge of collision cross sections is there- 
fore essential in understanding prospects for the success 
of either approach. 

We have previously considered electrostatic trapping of 
polar Il-state molecules from the point of view of stability 
with respect to collisions |l|| . The main bottleneck here 
is that electrostatic trapping requires the molecules to be 
in a weak-field-seeking state, in which case the molecules 
of necessity have an even lower-energy strong-field seek- 
ing state. Collisions involving the strong and anisotropic 
dipole-dipole interaction between molecules appear more 
than adequate to drive the molecules into these unfa- 
vorable states, leading to unacceptably high trap loss 
and heating. For this reason, it may be necessary to 
seek alternative methods that can confine the strong- 
field-seeking, absolute lowest-energy ground state of the 
molecule using time- varying electric fields 0, El El EE- 
A recent proposal for an electro dynamic trap is based on 
the microwave analogue of the familiar far-off-resonant 
optical dipole trap - only the microwave version doesn't 
have to be far off-resonance, making the trap very deep 

EE- 

Polar fermions may have an important advantage for 
electrostatic trapping, namely, low inelastic rates at cold 
temperatures. Kajita has discussed state- changing colli- 
sions of dipolar fermionic molecules |2l|, based on the 
well-known Wigner threshold laws for dipolar interac- 
tions. Namely, elastic scattering cross sections are essen- 
tially independent of collision energy E at low energies, 
but state-changing cross sections scale as E 1 ^ 2 . There- 
fore, at "sufficiently low" temperatures, elastic scatter- 
ing always wins, and evaporative cooling should be pos- 



sible. Using the Born approximation, Kajita concludes 
that this is the case for the molecules OCS and CH3CI, at 
reasonable experimental temperatures |2l| . This analysis 
may yet prove too optimistic, since the results include re- 
gions where the Born approximation may not be strictly 
applicable |2^] . Still, the idea is a sound one that deserves 
further exploration. 

A complete theoretical description of molecule- 
molecule scattering is complicated by the complexity of 
the short-range interaction between molecules. Indeed, 
for open shell molecules the potential energy surface is 
difficult to compute by ab intio methods, and remains 
inadequately known. It is therefore worthwhile to seek 
situations in which the influence of short-range physics 
is minimal. It appears that for weak- field seeking states 
the influence of the short-range potential is weak, o wing; 
to avoided crossings in the long-range interaction [lq . 
For collisions of identical fermionic molecules, the influ- 
ence of short-range physics may be even smaller, since 
only partial waves with I > 1 are present, and there is 
centrifugal repulsion in all scattering channels. 

A main aim of the present paper is thus to explore 
the suppression of inelastic collisions in fermionic Il-state 
molecules, using the OD radical as an example. This is 
an illustrative choice of molecule, since we have stud- 
ied its bosonic counterpart, OH, extensively in the past 
El HE 13 • ft 1S a l so a species at the center of current 
experimental interest [E EE • To this end we employ full 
close-coupling calculations to a model of the OD-OD in- 
teraction that includes only the dipolar part. We find, as 
we must, that the fermionic threshold laws ultimately fa- 
vor elastic over inelastic scattering at low temperatures. 
For OD, however, we find that the energy scales for this 
to happen remain quite low, on the order of microKelvins 
or below, so that the usefulness of this result to evapora- 
tive cooling remain questionable. 

On the bright side, the suppression of inelastic colli- 
sions does mean that a gas that is already cold will be sta- 
ble under collisions, even in an electrostatic trap. This is 
a similar conclusion to one we have drawn in the past for 
magnetostatic trappin g of spin-polarized paramagnetic 
(nonpolar) species [2a. l27j| . This is useful for cold colli- 
sions studies, since it is believed that collisions of weak- 
electric-field seekers are dominated by, and can be under- 
stood in terms of, purely long-range dipolar forces [23| . 
In particular, such collisions are predicted for bosons to 
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have long-range resonant states, termed field-linked res- 
onances, that may be useful in understanding cold colli- 
sions. A second goal of this paper is to verify that the 
fermionic OD molecules also possess these resonances. 



II. THRESHOLD LAWS IN THE BORN 
APPROXIMATION 

Threshold laws for various power-law long-range po- 
tentials have been written about extensively in the past 

IllllilMIMISIMHIllllillllMlllliilil^ in 

this section we summarize the main results relevant to 
the energy dependence of cross sections, using the first 
Born approximation to make the math transparent. Sim- 
ilar arguments are presented in Refs. I2H]. 

A first point to be considered is why the Born approxi- 
mation should be of any use at all, since it is ordinarily as- 
sociated with collisions of "fast" particles. Strictly speak- 
ing, however, the Born approximation is valid when the 
potential responsible for scattering is suitably "weak," 
meaning that the true scattering wave function is well- 
approximated by the unperturbed wave function. For 
dipolar scattering, the argument is as follows. Consider 
elastic scattering in a single- channel whose long-range 
potential varies as \/R s . Then, in partial wave I, the 
clastic scattering phase shift will vary with wave number 
k as 
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where a and (3 are constants depending on details of the 
potential. 

Thus for a dipolar potential with s — 3, the second 
term in Eqn. Q is the dominant contribution to the 
phase shift for any partial wave I > 1, yielding Si ~ k. 
(Moreover, the I = contribution to a realistic dipole- 
dipole interaction rigorously vanishes by symmetry.) It 
can be shown (for example, using the JWKB approxima- 
tion |4jj) that the second contribution in Eqn. (QJ arises 
from purely long-range physics, i.e., for intermolccular 
separations outside the centrifugal barrier imposed by the 
partial wave. As the collision energy approaches thresh- 
old, this distance gets ever larger, and the influence of 
the 1/i? 3 perturbing potential gets ever weaker. Thus, 
near threshold, the wave function is well-approximated 
by unperturbed spherical Bcsscl functions in each partial 
wave, and the Born approximation can be used. 

We adopt this view in the multichannel case. Because 
the dipole-dipole interaction is anisotropic, different par- 
tial waves are coupled together. Nevertheless, the diag- 
onal pieces of the Hamiltonian matrix have the general 
form 
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where R is the distance between molecules, \x is their 
reduced mass, and the effective C3 coefficient depends, in 



general, on the channel as well as on the degree of electric 
field polarization (see Ref . [l(| ) . When C| ff is negative, 
the potential J5J presents a finite barrier of height 
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For energies E considerably less than Eb, scattering only 
occurs from outside the barrier (barring resonances |IJ), 
thus setting an energy scale for the utility of the Born 
approximation. To make an estimate of this energy scale, 
consider the strong-field limit, where polarized molecules 
have C3 ~ d 2 , the square of the dipole moment. For OD, 
this sets the relevant p-wave centrifugal barrier height at 
■~ 10 nK. At higher energies, the incoming wave spills 
over the barrier, samples smaller-i? interactions, and is 
no longer well-described as a plane wave. 

Assuming that the Born approximation holds, we pro- 
ceed as follows. The partial scattering cross section for 
a collision process entering on channel i and exiting on 
channel / is given in terms of the transition matrix T by 



*if = ^\(i\T\f)\ 2 , 



(4) 



where the channel indices i and / include partial wave 
contributions U and If, which need not be the same. In 
the first Born approximation, the T matrix elements are 
given by the matrix elements of the potential (Chap. 7 
of Ref. |42j , where we have re-inserted the dimensionful 
factors) 



;//,/:= 2(^)(fc 2 fc/) 1/2 



(5) 



^ R^dRj h (hR)^^ Jlf (k f R). 

Here C^{lilf) represents the appropriate off-diagonal cou- 
pling matrix element, which, again, depends on field. 

For elastic scattering, where the initial and final wave 
numbers are equal, kf = fcj, we can rewrite Eqn. J5j in 
terms of the dimcnsionless variable x = kiR, 



(i\T\f) = 
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dx 
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The integral in Eqn. © converges whenever li + lt > 0, 
and is moreover independent of k%. Therefore, for any 
elastic scattering process by dipolar forces that changes 
I by at most 2 units, T ~ ki at low energies, and by 
Eqn. (@J, the associated cross section is independent of 
collision energy. In particular, the elastic scattering cross 
section of identical fcrmions does not vanish, if they in- 
teract via dipolar forces. 

For completeness, we give the value of the inte- 
gral. This is found by substituting ordinary Bessel 
functions for the spherical Bessel functions, j n (x) = 
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\J ir /2xJ n+ i/2(x), and using standard formulas for inte- 
grals [13]: 
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For an exothermic process, with A:/ > fcj, a similar 
argument yields for the transition amplitude 
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At 

dRJ li+1/2 {kiR)J lf+1/2 (k f R)R- 2 . 
This integral, too, can be done as long as U + If > [43| : 

/•oo 

/ dRJ li+1/2 (kiR)J lf+1/2 (k f R)R- 2 (9) 

^+l/2 r(ii + k) 
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where F stands for a hypergeometric function. 

Near threshold in an exothermic process, we have 
kf 3> fcj. In this case the leading order term of the 
hypergoemetric function F is a constant, and the only 
remaining dependence of @ on fej is in its prefactor. 
Thus T ~ fc' I+1 , and cr ~ fef i_1 - E 1 '- 1 / 2 . When the 
incident partial wave is = 0, as would be the case for 
identical bosons, the inelastic scattering cross section di- 
verges at threshold. For any higher partial wave, say the 
li = 1 partial wave that dominates scattering of identical 
fermions, the inelastic cross section instead vanishes in 
the threshold limit. 
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FIG. 1: Stark effect for the 2 H-3/2\f, rrif, parity) ground states 
of OD. Shown are state of f parity (a) and e parity (b). 



III. COLLISION CROSS SECTIONS FOR OD 

The OD radical differs from OH in two significant ways, 
for our present purposes: first, its lambda-doubling con- 
stant is somewhat smaller |44[ . Second, its hypcrfine 
structure depends on the nuclear spin of deuterium be- 
ing 1 instead of 1/2 for hydrogen, meaning that total spin 
states / = 1/2,3/2, and 5/2 are possible in the 2 n 3 / 2 
electronic ground state of OD (as in our OH work, we 
consider exclusively in the electronic ground state, and 
neglect excited vibrational and rotational levels). Figure 
1 presents the Stark effect for OD, which can be com- 
pared to the similar figure for OH, Fig. (1) of Ref.[l6|. 
Note that, due to the smaller lambda-doublet in OD, this 
radical enters the linear Stark regime at applied electric 
fields of ~ 200 V/cm, as opposed to - 1000 V/cm in OH. 



We consider collisions of the highest-energy weak- 
field seeking state in Fig.l, with quantum numbers 
|/, rrif, parity) = |5/2, 5/2, /). The details of our scatter- 
ing theory have been presented elsewhere, for OH (lt| . 
The main difference in handling OD is to incorporate 
Fermi exchange symmetry, which amounts to changing 
plus signs to minus signs in Eqn. (17) of Ref. 16]. Oth- 
erwise, we treat the scattering in the same way, by includ- 
ing only the Stark and dipole-dipole interactions, along 
with the hyperfine structure. 

Fig. 2 shows the main scattering results, as collision 
cross sections versus collision energy at different electric 
fields. Three different applied electric fields are indicated 
by color coding: £ — (black), £ — 100 V/cm (red), 
and £ — 1000 V/cm (blue). In each case, the solid line 
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FIG. 2: Elastic (solid lines) and total inelastic(dashed) 
cross sections for several values of electric field: £ — 
0, 100, WOOV/cm. 



denotes elastic scattering, while the dashed line repre- 
sents total inelastic scattering to channels where one or 
both molecules loses internal energy. In zero field, the 
molecules are completely unpolarized, and the dipole- 
dipole interaction vanishes. Thus the cross sections obey 
the familiar Wigner threshold laws for short-ranged in- 
teractions between fermions: the elastic cross section 
a e i oc E 2 , whereas the (exothermic) inelastic cross section 
dinei oc E 1 / 2 jljj. Thus in zero field elastic scattering is 
actually less likely than inelastic scattering at lower en- 
ergies (below about in this example). Above this 
energy elastic scattering appears somewhat more favor- 
able than inelastic scattering, at least until several mK, 
where both cross sections start to hit the unitarity limit. 

Turning on the electric field partially polarizes the 
molecules, so that the dipole-dipole interaction is "acti- 
vated." Then the dipole-dipole threshold laws take effect: 
a e i cx const., whereas we still have Ui ne i oc E 1 / 2 . Fig. 2 
illustrates where this threshold behavior kicks in for dif- 
ferent electric field values. Notice that the higher the 
electric field, the lower is the energy where the thresh- 
old behavior is attained. This is because the effective 
C3 coefficient that determines the barrier height Ef, is 
an increa sing function of electric field, at least until it 
saturates |l6|. 

On the other hand, a Fermi gas of molecules that is 
already cold will enjoy the benefits of Wigner-law sup- 
pression of inelastic collisions. Suppose a quantum de- 
generate gas of OD could be produced at nK tempera- 
tures, as is the case for current experiments in 40 K and 
6 Li. Then Fig. 2 suggests that a small bias field of ~ 100 
V/cm reduces inelastic cross sections to an acceptable 
level of ~ 2 x 10~ 14 cm 2 , corresponding to a rate con- 
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FIG. 3: Comparison of cross sections for OD (red) and OH 
(black) molecules. Solid and dashed lines refer to elastic and 
inelastic cross sections, respectively. 

stant ~ 10 -16 cm 3 /sec. 

To emphasize the difference between bosons and 
fermions, we reproduce the £ = 100 V/cm cross sec- 
tions in Fig. 3 (red) along with the corresponding cross 
sections for OH in the same field (black). It is clear that 
in both cases elastic scattering (solid lines) is quite simi- 
lar, whereas the behavior of inelastic scattering is utterly 
different at low energies for the two species. Equally im- 
portantly, at collision energies above about 1 mK, all the 
cross sections have the same general behavior. This is a 
manifestation of the strength of the dipolar interactions, 
and the fact that in this energy range all processes are 
essentially unitarity-limited. 



IV. ON THE QUESTION OF FIELD-LINKED 
RESONANCES 

Finally, we comment on the occurrence of field-linked 
(FL) resonance states in this system. Fig. 4a shows the 
elastic and inelastic cross sections versus electric field, at 
a fixed collision energy of 1 /xK. This figure exhibits the 
characteristic peaks indicative of field-linked resonances; 
compare Fig. (2) of Ref. [3. To converge these results 
at higher field demands an increasing number of par- 
tial waves. Fig. 4b illustrates the convergence of the 
resonant scattering cross section for various numbers of 
partial waves included. For partial waves L = 1,3,5,7, 
the cross section is well-converged up to several hundred 
V/cm. This is sufficient to compute the first two reso- 
nance states, which are the only well-resolved ones any- 
way. 

As discussed in Ref. 0] and elaborated on in Ref. 
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for s- waves, and makes an effect only at second order 
llfl . For fermions in identical spin states, however, the 
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FIG. 4: a) Elastic (solid) and inelastic (dashed) cross sections 
for OD scattering as a function of applied electric field. The 
Collision energy is E — lfiK. b) Convergence of elastic cross 
section upon increasing the number of partial waves included 
in the calculation. 



[24J , the FL resonances live in adiabatic potential energy 
surfaces generated by avoided crossings in the long-range 
dipole-dipole interaction. The surfaces have a somewhat 
different character for fermions than for bosons, however. 
Fig. 5 shows a set of adiabatic curves for a single partial 
wave L = 1; this is the analogue of Fig. (1) in Ref. |23j . 

For identical bosons, the long-range attractive part of 
the relevant curve is predominantly s-wave in character, 
hence has a 1/R 4 behavior in an electric field. This re- 
flects the fact that the direct dipolar interaction vanishes 



FIG. 5: A set of avoided crossings that generate FL states for 
OD. 



attractive part involves the p-wave interaction, where the 
dipole is already nonzero, so that the long-range inter- 
action scales as I/i? 3 . The net effect is that the inner 
turning point of the s-wave FL states approaches smaller 
R as the field is increased for bosons, but that this inner 
turning point is relatively fixed for fermions. (The outer 
turning point is set by the energy of the resonant state 
relative to the threshold, and is thus arbitrarily large.) 

These resonant states, if sufficiently stable, may form 
a novel kind of pair of fermions, which may ultimately 
lead to an exotic Fermi superfluid state. Unfortunately, 
as seen in Fig. 4, these resonances are quite readily sus- 
ceptible to predissociation, indicated by the large inelas- 
tic cross sections near resonance. In this they resemble 
their bosonic counterparts. However, stabilization of cold 
dipolar gases using magnetic fields has been recently dis- 
cussed 45] . It is yet conceivable that these resonances 
could be tamed long enough to put them to use. 

In summary, we have computed scattering cross sec- 
tions for cold collisions of the fermionic free radical OD, 
as functions of both collision energy and electric field. 
We find that, similar to the case of bosonic OH, these 
molecules are unlikely to be stable against collisions in 
traps warmer than about 10 /iK. Unlike OH, however, 
they will be collisionally stable at lower temperatures, 
owing to the unique Wigner threshold laws for fermionic 
polar particles, in such a gas. 

This work has been supported by the NSF and by a 
grant from the W. M. Keck Foundation. We acknowledge 
useful interactions with M. Kajita and D. DeMille. 
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